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„  .  SUMMARY 

M*  know-ln  lSle^ly^a  positive  and  a  negative  value  of  a 
function  at  two  given  points  respectively.  The  function  Is 
continuous  and  convex  and  Is  otherwise  unknown  but  computable. 

v  Sj- 

Starting  with  this  Information,  ws  dosorlbe  a  procedure^ for 
locating  its  unique  root  (on  the  starting  Interval)  within 
an  Interval  of  minimum  guaranteeable  length  In  n  steps,  where 
a  step  consists  of  calculating  the  value  of  the  function  at 
any  point  we  choose .  The  pertinent  functional  equation  is 
derived  and  curves  of  the  objective  function  are  plotted  for 
n  -  1,  2,  3,  4  from  data  obtained  from  the  Johnnlac,  RAND’s 
Princeton— type  high-speed  digital  computer.  (  )  _ 
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1 .  INTRODUCTION 

In  the  past,  methods  have  been  proposed  and  used  for  the 
computation  of  real  roots  of  equations  in  a  single  variable. 
Notable  among  these  is  Newton’s  Tangent  Method.  However,  this 
procedure  is  quite  costly  in  many  instances  wherein  the  function 
whose  zero  is  desired  is  inordinately  complicated,  and  Indeed, 
it  is  usually  the  case,  polynomials  and  other  simple  functions 
excepted,  that  derivatives  are  more  complicated  to  evaluate 
than  the  functions  themselves.  This  last  circumstance  con¬ 
stitutes  the  principal  motivation  of  later  methods  which  are 
based  solely  upon  evaluations  of  the  functions  themselves.  We 
shall  not  describe  any  of  these  in  detail  here,  however,  but 
shall  content  ourselves  with  a  single  reference  [l] ,  and  the 
remark  that  our  method  falls  in  this  latter  category.  To  our 
knowledge,  none  of  these  later  methods  are  sequential  mlnimax 
in  character  except  for  the  bisection  procedure  and  this  in  a 
restricted  sense. 

Let  us  now  mention  that  due  to  personal  preference,  or 
otherwise  established  physical  laws,  the  functions  encountered 
in  the  past  and  presumably  to  be  encountered  in  the  future,  had 
a  strong  tendency  to  be  analytic  in  character.  This  happy 
circumstance  coupled  with  the  rare  occurrence  of  the  zero 
sought  being  a  point  of  inflection  of  the  graph  of  the  function 
Justifies  the  assertion  that  once  a  root  is  located  on  a  suf¬ 
ficiently  srall  Interval  by  a  positive  and  negative  reading, 
the  function  will  be  inclined  to  be  convex  or  concave  throughout 


the  interval. 
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Thua,  for  example,  in  [l]  the  equation  f(t)  ■  t*-6ta 
+9t+6-6tae”t-l8te“’t-3te“2t-6ta  »  0  is  studied.  f(2)  and  f(3) 
were  found  to  be  of  oppoaite  sign  and  thua  a  root  waa  located 
on  (2,3).  By  repeated  application  of  the  "bisection"  technique 
the  root  waa  found  to  be  approximately  2.7184.  However,  to 
obtain  4  decimal  accuracy  with  certainty  ualng  the  blaectlon 
technique  would  almoat  alwaya  require  (Btartlng  on  an  Interval 
of  unit  length  aa  above)  14  evaluations  of  the  function  even 
with  the  Information  that  the  function  la  monotone.  On  the 
other  hand,  it  can  be  shown  by  elementary  arguments  that  the 
above  function  f  la  convex  on  [2,3]  •  Making  use  of  this  in¬ 
formation,  together  with  the  knowledge  of  f(2)  and  f(3)  and 
succeeding  functional  values  (the  blaectlon  procedure  makes  use 
only  of  the  knowledge  (x^,  agn  f(x^))  one  should  expect  to 
obtain  far  mor*.  rapid  convergence  using  a  sequential  minimax 
procedure  which  takes  this  information  into  account.  This 
phenomenon  will  be  exemplified  in  paragraph  3. 

2.  DESCRIPTION  OF  THE  PROCSDPKS 

In  what  follows  we  shall  describe  a  numerical  procedure 
for  solving  the  following  problem: 

"We  know  initially  a  positive  and  a  negative  value  of  a 
function  at  two  given  points  respectively.  The  function  is 
continuous  and  convex  and  is  otherwise  unknown  but  computable. 
Given  an  Integer  n  >  0,  how  do  we  proceed  to  locate  its  root 
on  this  initial  interval  within  an  interval  of  minimum  length 
in  n  steps  where  a  step  consists  in  calculating  the  value  of 
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the  function  at  any  point  we  choose?"  As  It  stands,  the  ques¬ 
tion  is  meaningless,  i.e.,  has  no  definite  answer,  until  we 
specify,  for  example,  that  our  procedure  be  sequential  rain— 
irr.ax,  i.e.,  to  say  at  each  step  of  the  procedure  we  assume  that 
the  worse  possible  situation  might  occur  from  that  point  on  in 
the  light  of  our  present  information  about  the  function  and 
proceed  to  evaluate  the  function  at  a  point  which  hedges  against 
all  contingencies  so  as  to  guarantee  a  fixed  interval  at  that 
stage . 

An  alternate  question  would  envolve  the  final  interval  as 
specified  in  length  and  the  requirement  that  the  number  of  steps 
be  minimized.  Without  going  into  any  detail,  we  mention  that 
the  procedure  which  answers  the  first  question  need  be  somewhat 
modified  to  answer  the  second. 

Plnally  we  remark  that  if  the  function  is  concave  rather 
than  convex,  we  merely  have  to  replace  the  function  by  Its 
negative,  i.e.,  f(x)  — >  —  f(x)  for  the  procedure  to  be  applicable. 
Also,  if  the  signe  of  the  initial  readings  are  reversed,  i.e., 
a  <  b  and  f(a)  <  0,  f(b)  >  0,  where  (a,b)  is  our  starting  inter¬ 
val,  we  need  only  replace  f(x)  by  f(— x)  to  render  the  procedure 
applicable . 


We  now  describe  the  procedure  cycle: 

Suppose  we  are  in  the  situation  in  which  we  know  f(a)  -  Yft>  0, 
r(b)  -  ~Y^,  Y^  >  0,  where  a  <  b,  and  the  root  >  S  and  we  have 
n  more  readings  to  make . 


W  -  a  ♦  (b-a) 
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If  n  ■  0,  the  computation  ceases  and  the  values  S,W  are 
recorded . 


If  n  >  0,  calculate  the  value  of  f  at  x  ■  S  +  (W-S)*/°  (w-)  * 

xa 

Y_ 

If  f(x)  ■  Y  >  0,  set  a*  •  x,  b’  •  b  and  S'  ■  x+(xr-e)  y  _y-  . 

x  xa  Xx 

If,  however,  f(x)  -  -Y„  <  0,  set  a*  -  a,  bj  ■  x,  and  if 

x  Y 

Y  >  Yh  set  S'  -  S,  otherwise  set  S'  •  max(S,x-(b-x)  « — £-  ). 

*  D  xb  xx 

Finally,  set  n‘  -  n-1 . 

Then  we  are  in  the  situation  in  which  we  know  f(a')  •  Yftl>  0# 

f(b')  •  -Ybi#Yb,>  0,  where  a'  <  b'  and  the  root  >  S'  and  we 

have  n'  more  readings  to  make.  This  completes  the  cycle.  (As 

the  problem  is  stated,  S  •  a  initially.) 

In  the  next  section  we  shall  illustrate  how  this  procedure 

works  on  a  particular  example. 

Remark  1 .  The  foregoing  procedure  is  an  approximation  to 

the  actual  minimax  procedure.  The  theoretically  correct  pro— 

/Yb\ 

cedure  would  only  involve  replacing  the  expression  Pn \y~) 

ft 

Y 

the  formula  for  x  above  J^)  *  wilere  ^n(S,Y)  is  defined 

81 

in  paragraph  4.  Since  the  objective  function  is  relatively  in¬ 
sensitive  to  S  in  our  choice  of  /°n(S,Y)  in  the  vicinity  of  the 
minmax,  we  ftsl  that  the  approximation 

/*n(0,Y)  -  /®n(S,Y) 

is  Justified,  and  define  ^n(Y)  ■/3n(0,Y).  This  approximation 
renders  the  procedure  more  adaptable  to  machine  computation. 


paper . 


#Qraphs  of  p  for  n  -  1,2,3#**  are  included  at  the  end  of  this 


Remark  2 .  Since  the  average  digital  computer  has  difficulty 


in  reading  graphs,  in  order  to  program  the  procedure,  one  would 
first  find  suitable  algebraic  approximations  to  the  graphs  of 
/°n  after  the  fashion  of  C.  Hastings,  Jr.  [}) . 

A  NUMERICAL  APPLICATION— COMPARING  THE  BISECTION  TBCHN1QPB 

Suppose  we  are  desirous  of  bracketing  in  on  the  zero  of  a 

certain  complicated  function  f  defined  over  the  interval  (0,1), 

We  know  that  f  is  continuous  and  convex  and  that  in  fact,  f(0)  *1, 

f(l)  -  -1.  However,  since  the  function  requires  one  hour  of 

machine  time  to  evaluate  a  single  point,  we  are  ignorant  of  the 

fact  that  to  all  intents  and  purposes  it  is  given  by  the  rel— 

atively  imocuous  expression  f(x)  -  max(-l,(x-  *><? 

can  afford  to  make  three  evaluations  of  the  function  fn-2). 

Y0 

Upon  referring  to  the  graph  of  R»(0,Y)#  with  Y  w  — ^  ■  1  we  see 

xi 

that  we  can  guarantee  locating  the  root  on  an  interval  of  length 
.01  times  the  original  interval,  i.e.,  on  an  interval  of  length 
.01.  However,  since  the  graphs  represent  the  worst  that  can 
happen  to  us,  we  expect  to  do  much  better,  and  indeed  this  turns 
out  to  be  the  case. 

We  proceed  to  calculate: 

Cycle  Ho  .  1 . 

a  -  0,  b  -  1,  Ya-  1,  Yb  -  1,  S  ■  0,  n  ■  3,  whence  by  our 
formula,  W  -  and  we  have  located  the  root  on  (0,.^).  Next, 

- i - 

At  the  end  of  this  paper. 


P-935 

9-n^ 

x  -  0  +  (.5-0)  /»(1)  *  .5(.l48)  -  .074,  and  we  find  that 
f(.074)  -  .76839  >  0;  so  a'  -  .074,  b*  -  1  and  •  •  -  .074 
*  *  *51950.  Finally,  n’  -  2. 

Cycle  No.  2.  (dropping  primes  on  the  new  variables) 
a  -  .074,  b  -  1,  Ya  -  .76879,  Yb  -  1,  S  -  .31950,  n  -  2, 
whence  by  our  formula,  W  ■  .074  +  •  .47636,  and 

we  have  located  the  root  on  (.31950,  .47636).  Next, 

X  -  .31950  +  .15686  .  .31950  +  (.15686)  (.198)  -  .35066, 

and  we  find  that  f( .35066)  -  -.04795  <  0,  ao  a’  -  .074,  b'  -  .35066 
and  since  Yx  <  Tfa,  S'  -  max  (.31950,  .35066  -  (1-^066|^047^)  ) 

-  .31950.  Plnally,  n'  -  1. 

Cycle  Wo.  3. 

a  -  .074,  b  -  .35066,  Yfi  -  .76839,  Yfe  -  .04795,  S  -  .31960, 
n  -  1,  whence  by  our  formula,  W  -  .074  ♦  ^  .33441 

and  we  have  located  the  root  on  (.31950,  .33441).  Next, 
x  -  .31950  +  ( .33441-. 31950)  {^)  -  .31950  +  (.01491)  ^(.624) 

-  .32440,  and  find  that  f ( .32440)  -  .02535  >  0,  ao  a'  -  .32440, 
b'  -  .35066,  and  S'  -  .32440  +  fyUffyrv ^ ^25221  -  *33294. 

Plnally,  n*  *0. 

Cycle  No.  '4 . 

a  -  .32440,  b  -  .35066,  Y&  -  .02535,  Yb  -  .04795,  S  -  .33294, 
n  •  0,  whence  by  our  formula,  W  -  .32440  ♦  ^ 

■  .33348,  and  we  have  located  the  root  on  ( .33294 ,  .33348)  •  n  •  0, 


so  the  computation  ceases  and  the  Interval  (.33294,  .33348)  is 
recorded . 


Discussion 
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Slnce  the  bisection  procedure  does  not  take  c< gnizance  of 
the  convexity  of  the  function,  one  would  obtain  initially  without 
any  evaluation  of  f  that  the  root  Ilea  on  (0,1);  with  one  evalua¬ 
tion,  (0,.5)»  with  two  evaluatlone,  (.25*. 5);  and  finally,  with 
three  evaluations,  (.25*. 375).  The  lengtha  of  these  bracketing 
intervals  are,  then,  1,  .5,  .25,  .125  respectively.  Comparing 
these  with  those  obtained  by  the  procedure  above,  namely  .5, 

.15686,  .01491,  .00054,  we  see  that  our  procedure  does  have  a 
slight  edge  in  this  instance.  In  fact,  in  any  instance  of  a 
convex  function  with  the  same  starting  values  as  our  example, 
we  can  guarantee  in  three  evaluations,  a  bracketing  interval 
of  length  ^  .01,  as  was  pointed  out  earlier.  This  value  compares 
favorably  with  .125. 

4.  DERIVATION  OF  THE  FUNCTIONAL  EQUATION 

Given  a  continuous  convex  function  f  and  two  of  its  values, 
one  positive  and  the  other  negative,  (Bay  f(xt)  >  0  and  f(x*)  <  0 
with  X|  <  Xu)  how  should  one  search  for  the  zero  of  the  function 
which  lies  on  the  known  interval  (xi,x*)?  Utilizing  the  principle 
of  optimality  of  the  theory  of  Dynamic  Programming  [2]  ,  we  formulate 
the  problem  in  terms  of  minimizing  the  maximum  length  of  interval 
on  which  we  can  deduce  that  the  zero  is  located  in  n  readings 
taken  sequentially. 

First,  by  reduction  to  scale,  we  can  always  consider  the 


following  diagram: 
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1 

k 


Fig.  1 


where  f (0)  ■  1,  f{l)  -  -Y,  Y  >  0,  W  •  yiy.  Since  f  ia  convex, 
the  eero  must  lie  on  (0,W]  .  If  we  have  Just  one  more  reading 
to  make  (n-l),  then  we  choose*  x  on  (0,W)  and  calculate  f(x). 
It  can  be  shown  by  simple  dominance  arguments  that  no  reading 
of  f  need  ever  be  taken  outside  any  interval  on  which  the  zero 
has  been  located.  Having  chosen  x  and  calculated  f(x)  we  en¬ 
counter  exactly  one  of  the  following  oases  (barring  f(x)  •  0, 
of  course,  the  best  possible  case). 

Cast  1 

If  f(x)  -  v  >  0,  then  by  drawing  straight  lines  Joining 
- 1 - 

The  particular  optimal  choice  will  be  derived  in  what 
follows . 
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known  points  on  the  graph  of  f,  we  have  the  picture  below,  with 
the  root  located  on  [S,W]  as  implied  by  the  convexity  of  f. 


1 


1 


Y 


Fig.  2 


Case  2 

If  f(x)  •  — v '  <  0,  the  picture  is 
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Pig.  5 


with  the  root  located  on  (max  (0,t'),W")  as  Indicated. 

Let  S"  ■  max  (0,S’).  If  x  la  the  final  point  at  which 
f(x)  ia  to  be  determined  (l.e.  n»l)  then  it  ia  chosen  so  as 
to  minimize  the  maximum  possible  value  of  max  (W'-S,  W,1-Sr1) 
consistent  with  our  choice. 

Before  going  Into  the  algebraic  details  of  the  solution 
for  n«l,  let  us  consider  the  general  n  stage  process  in  which 
we  have  several  more  readings  to  make.  In  either  of  the  two 
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cases  diagrammed  above,  all  the  essential  data  can  be  described 
by  a  basic  triangle  determined  by  a  two  parameter  system  as 
follows  . 

In  Pig..  2,  we  can  conclude  that  the  graph  of  f  lies  above 
the  line  segment  vS  and  below  the  segment  vW’Y  (where  the  letters 
stand  for  both  points  and  values  in  the  obvious  manner).  If 
we  now  draw  the  line  segment  SY,  it  may  or  may  not  be  crossed 
by  the  graph  of  f;  but  if  it  is  crossed  at  some  point  P  say, 
we  can  replace  that  portion  of  the  curve  lying  below  PY  by  PY 
with  impunity.  This  device  does  not  violate  the  convexity 
condition  nor  does  it  exclude  the  worst  possible  case.  This 
can  be  shown  by  simple  dominance  argument*  based  on  our  knowledge 
of  the  function  at  any  particular  stage. 

The  essential  data  can  thus  be  described  by  the  triangle 
vSY.  Since  a  vertical  reduction  in  scale  leaves  the  problem 
invariant,  and  a  horizontal  reduction,  relatively  invariant,  the 
triangle  vSY  can  be  replaced  by  a  triangle  of  standard  form 
described  pictorial  ly  by  two  parameters  Y,S  (say)  as  follows. 
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Fig.  4 

Similarly,  in  the  second  case  (Fig.  5)#  the  graph  of  f 
lies  above  S'v'  and  below  lW"v'.  Draw  the  line  segment  from  1 
to  S'  and  replace  any  portion  of  the  graph  of  f  lying  below 
this  line  by  the  line  itself  from  1  to  the  point  of  crossing. 
This  does  not  effect  the  choice  of  subsequent  x's,  since  they 
will  all  be  chosen  on  minimal  bracketing  intervals.  Again, 
this  can  be  shown  by  simple  dominance  arguments.  Thus,  in  the 
second  case  also  by  a  suitable  reduction  to  scale  we  are  led 
to  another  representation  of  Fig.  4. 

Now  define 

R^SjY)  •  the  minimum  length  of  interval  on  which  we  can 
guarantee  locating  the  zero  in  [0,l]  of  any  convex  function  f, 
given  that  f(0)  ■  1,  f ( 1 )  ■  -Y  <  0,  the  root  is  >  S  and  we  have 
n  readings  to  make.  If  n  -  0,  we  have  clearly. 
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rv(stY)  - 


1 

T+Y 


-  S 


Next,  using  the  principle  of  optimality,  and  taking  into  account 
the  scale  factors,  we  obtain  for  n  >  0  the  following  recurrence 
relation : 


Rn(S,Y)  -  min  ,  max 


XVl  (%£rj  *v'j 


max  (l—x)  Py,  1  (  r^7  *  T“7 

0<v^-x(UY)  n~l  V  1_x  ^ 


with  the  upper  and  lower  expressions  after  the  brace  correspond¬ 
ing  to  the  second  and  first  cases,  respectively.  The  scale 
factors  are  obtained  in  a  completely  elementary  manner  by  means 
of  similar  triangles. 

The  ranges  of  the  variables  S,Y  above  are  given  by  Y  2  0, 

0  £  S  £  .  To  render  the  expressions  more  amenable  to 

computation  on  the  Johnniac,  we  made  the  following  substitutions  1 


W  -  yiy  ,  <|>n(S,W)  •  1^(3, Y),  whence  (S ,  Y )  -  in(S»y^y)» 

An  additional  modification  of  the  quant  led  variables  v,  v* 
finally  reduces  the  system  to 


*0(S,W)  -  '^S 


4  { S  ^  ■  min  max 

S<x<W 


(1-x)  max  4  . 

x^t’<W  rwi 


t  '-x 
T=x” 


W(t*-x 
VOJx  ' 


where  0  <;  S  <;  W  ^  1 . 
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The  functions  4>n(S,W)  were  then  computed  for  n  •  1,2, 3, 4 
by  means  of  a  discrete  approximation  using  various  grid  sizes 
and  linear  interpolation.  We  are  indebted  to  S.  Dreyfus  for 
programming  these  computations  for  the  Johnniac . 

The  minimizing  x’s  were  of  course  recorded  and  these  form 
the  basis  of  our  optimal  policy,  i.e.  to  say  x#-  x*(S,W)  is 
the  point  at  which  we  evaluate  our  unknown  function  f  given 
that  the  root  lies  on  (S,W)  in  our  basic  triangle  and  there 
are  n  evaluations  to  be  made.  The  point  x*  is,  of  course,  itself 
measured  on  the  basic  triangle.  To  take  c?re  of  the  general 
situation,  we  must  of  course  relate  our  readings  to  the  original 
scale.  If  we  let/^(8,Y)  denote  the  fraction  of  the  distance 
between  S  and  W  occupied  by  x*,  we  readily  obtain 


(s,w)— s 

— W=S - 


where  W  «  yyy  •  It  is  now  a  relatively  easy  matter  to  relate 
x*  to  the  original  scale  and  thus  to  obtain  the  procedure  cycle 
outlined  in  paragraph  2. 

Qraphs  of  the  functions  Rn(0,Y)  and  /?n(0,Y)  for  n  -  1,2,), 4 
are  Included  at  the  end  of  this  paper. 


S.  REMARKS — THE  SPECIAL  CASE  n  -  1 

We  shall  initiate  this  section  with  a  few  remarks  Intended 
primarily  to  validate  certain  assumptions  made,  tacitly  or  other¬ 
wise,  in  the  derivation  of  the  functional  equation.  We  shall 
close  this  section  with  a  brief  treatment  of  the  special  case 
n  »  1,  S  •  0,  which  provided  us  with  an  excellent  check  on  the 
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validity  of  the  data  obtained  from  the  Johnnlac .  The  remark® 
are  as  follows. 

1.  Suppose  we  are  at  a  certain  stage  In  an  optimal  sequential 
mlnlmax  search.  We  have  several  values  of  f  computed  and  we  are 
ready  to  choose  our  next  point  of  evaluation.  Let  a,b  denote 

the  closest  evaluation  points  on  the  left  and  right  respectively 
of  the  minimal  Interval  (S#W)  on  which  the  root  is  known  to  lie. 
Then  In  an  optimal  procedure,  we  choose  x  on  (a,b)  .  To  see 
this,  we  need  only  state  that  since  the  unknown  function  f  may 
indeed  be  piecewise  linear  (a  possible  contingency)  outside 
(a,b)  is  Is  easy  to  see  that  any  such  subsequent  reading  would 
in  such  a  contingency  afford  us  no  information  regarding  the 
character  of  f  within  (a,b)  pertinent  to  the  location  of  Its 
zero,  not  already  Implied  by  the  quantities  a  ,b,f (a ) , f (b) ,  S  and 
W,  or  indeed  by  any  of  our  previous  readings  for  that  matter. 
Similar  but  slightly  more  involved  arguments  can  be  given  to  show 
that  the  next  reading  ahould  be  taken  on  the  interval  (S,W). 

2.  In  the  treatment  of  Case  2,  it  was  tacitly  assumed  that 
Y  <  Y,  i.e.,  to  say  the  SY  line  has  a  negative  slope  as  shown 

A 

in  the  figure.  Again,  it  can  be  shown  by  dominance  arguments 
that  the  worst  situation  occurs  when  f  is  monotone,  and  Indeed 
when  the  graph  of  f  lies  abc  e  the  line  SY.  It  is  this  condition 
which  determines  the  limits  of  the  quantified  variable  v  in  the 
upper  line  of  the  functional  equation  for  Rn .  It  is  precisely 
such  dominance  considerations  as  these  which,  though  enabling 
us  to  express  the  functional  equation  in  a  relatively  simple  form 
and  to  obtain  an  optimal  "policy”  therefrom  via  its  recursive 
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computation  ,  do  not  enable  us  to  obtain  an  optimal  "procedure” 
directly  from  the  equation.  However,  all  the  happy  contingencies, 
together  with  the  information  about  the  function  f  they  afford 
are  taken  care  of  in  the  procedure  outlined  in  paragraph  2.  It 
is  intuitively  clear  that  the  resulting  procedure  is  sequential 
minimax . 

We  now  conjecture  that  the  functional  equation  for 
may  be  validly  simplified  by  assuming  that  the  maximum  in  the 
upper  line  of  the  equation  is  always  taken  on  at  the  upper  end- 
point  of  the  range  of  v.  This  is  true  for  example  if  n-1 .  However, 
the  authors  were  unable  to  prove  the  assertion  in  general  and 
thus  the  transformed  equation  for  was  subsequently  submitted 
to  the  Johnniac  in  its  present  form.  Suffice  it  to  say  that 
the  resulting  data  supported  the  conjecture. 

4.  R^SjY)  is  separately  decreasing  in  S  and  Y.  To  see 
this  for  the  first  variable,  for  example,  upon  recalling  the 
definition  of  ^(SjY),  we  need  only  observe  that,  other  things 
being  the  same,  the  information  that  the  root  is  >  S  Includes 
the  information  that  the  root  is  >  S'  if  S*  <  S.  On  the  basis 
of  the  additional  information,  then,  we  can  clearly  guarantee 
at  least  as  short  a  final  bracketing  interval  with  a  larger  S 
value  if  we  are  proceeding  optimally,  i.e.  R^SjY)  <  F^(S'#Y) 
if  S  >  S',  kn  analogous  argument  applies  for  the  second  variable. 

3.  We  conclude  our  remarks  with  a  conjecture  about  the 
asymptotic  behavior  of  ^(SjY)  as  Y  — >  ^  -  1 .  In  particular, 
we  conjecture  that  there  exists  a  positive  number/?*  such  that 
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/°  (O.Y)  J,  P*  as  Y  — >oo  .  This  la  easily  shown  to  be  true  for 
n  -  1 ,  but  the  proof  seems  to  be  difficult  for  larger  values 
of  n.  The  precise  value  for  n  ■  1  Is  known  and  Is  Included 
in  the  accompanying  graph  together  with  "estimates "  of  /°*  for 
n  -  2,3,4  . 

We  conclude  this  paper  with  a  brief  discussion  of  some 
results  obtained  for  the  case  n  *  1,  S  »  0.  We  shall  spare 
the  reader  the  elementary  albeit  complicated  algebraic  details 
involved  and  simply  state  that  upon  substituting  the  function 

Ro(5,Y)  -  y~y  ~  s 

in  the  right  member  of  our  functional  equation  and  setting  S  -  0 
in  the  result  we  obtain  a  relatively  simple  algebraic  minima* 
problem  for  the  determination  of  R|(0,Y). 

The  substitution  -  x(lV/)  then  yields,  upon  performing 
the  mlnimax  operation, the  following  equation  relating  the  optimal 
ratio  f'  and  the  variable  Y: 

+1  )Ya  ♦  2{2f*-f*~y>+\)X  ♦  {ZP-\  )*  -  0  , 


with  the  restrictions  Y  >  0,  0  <>  f>  <  1. 

By  some  inexplicable  quirn  of  fate  the  discriminant  of  the 
above  quadratic  In  Y  turns  out  to  be  precisely  8/^.  Thus,  we 
readily  obtain  the  following  rational  parameterization  of-  the 
(l^,Y)  curve: 


r-  ±(i-oa 

Y  -  4 1* 

1— 4 1— t4 


0  <  t  <  tQ 


J 
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To  determine  the  value  of  t^  we  note  that  the  limiting 
form  of  the  polynomial  equation  as  Y  — Xd  is  simply  the  co¬ 
efficient  of  the  leading  term  In  Y  set  equal  to  zero: 

2f*- bf*- -  0. 

This  equation  has  a  unique  root  on  the  interval  (0,1)  and  is 
given  by 


2 


This  is  the  asymptotic  value  mentioned  earlier  in  remark  5.  It 
follows  from  the  parametric  representation  of  !°  that  tQ  »  1  -j2?r 
to  .250.  As  a  check  on  this  we  note  that  this  value  of  t^  is 
the  smallest  positive  root  of  the  denominator  in  the  parametric 
expression  for  Y. 

The  graphs  of  Ri(0,Y),  /^(Y)  were  plotted  from  manual  com¬ 
putations  using  the  preceding  formulas  and  these  compared  quite 
favorably  with  the  results  from  the  Johnniac.  Unfortunately, 
however,  due  presumably  to  the  choice  of  insufficiently  small 
grid  sizes  imposed  by  the  Johnniac' s  limited  memory,  a  cumulative 
error  caused,  in  effect,  a  gradual  upward  creep  in  the  tails  of 
the  f>  curves .  These  were  smoothed  out  to  agree  with  theory  to 
obtain  the  included  graphs.  The  R  curves,  fortunately,  seemed 
quite  insensitive  to  choice  of  grid  si/,e. 

In  closing,  we  mention  that  further  checks  on  the  computa¬ 
tions  were  provided  by  the  easily  derived  relationships: 

+n(W,W)  -  0.  *n(S,l)  -  ip  , 
and  these  were  found  to  fit  the  Johnniac  data  exactly. 
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